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Abstract. Using a sample of 57 VLT FORS spectra in the redshift range 1.37 < z < 3.40 (selected mainly from 
the FORS Deep Field survey) and a comparison sample with 36 IUE spectra of local (z ~ 0) starburst galaxies 
we derive Civ and Siiv equivalent width values and estimate metallicities of starburst galaxies as a function 
of redshift. Assuming that a calibration of the C IV equivalent widths in terms of the metallicity based on the 
local sample of starburst galaxies is applicable to high-z objects, we find a significant increase of the average 
metallicities from about 0.16 Zq at the cosmic epoch corresponding to z ~ 3.2 to about 0.42 at z w 2.3. A 
significant further increase in metallicity during later epochs cannot be detected in our data. Compared to the 
local starburst galaxies our high-redshift objects tend to be overluminous for a fixed metallicity. Our observational 
results are in good agreement with published observational data by other authors and with theoretical predictions 
of the cosmic chemical evolution. 
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1. Introduction 

Among the unsolved questions of cosmology is the begin- 
ning and evolution of the star formation process at early 
cosmic epochs. From the presence of heavy nuclei in high-z 
quasars and galaxies it is clear that star formation started 
rather early (see e.g. Hamann and Ferland 1999, Dietrich 
et al. 1999). Moreover, galaxy counts and emission line 
studies of high-z galaxies indicate that the star formation 
rate (SFR) declined rapidly since at least the epoch cor- 
responding to redshifts of z ss 2 (see e.g. Madau et al. 
1996, Madau 2000). On the other hand intrinsic interstel- 
lar extinction (which may affect strongly the rest-frame 
UV of distant galaxies) and uncertainties concerning the 
relation between line emission and the star formation in- 
troduce considerable uncertainties into the derivation of 
the global star formation rate at early cosmic epochs from 
visual galaxy counts and emission line studies, while in- 



tend offprint requests to: D.Mehlert, 
e-mail: dmehlertOlsw . uni-heidelberg . de 

* Based on observations obtained with the FORS 
Instruments at the ESO VLT, Paranal, Chile on the course 
of the observing proposals 65.O-0049(A), 66.A-0547(A), 
68.A-0013(A), 68.A-0014(A) 



vestigations using IR and mm-wave data are hampered by 
small samples. 

An alternative approach to investigate the early star 
formation history of the universe is the evaluation of the 
chemical enrichment history of the universe, as star for- 
mation and the rapid evolution of massive stars results in 
a production of heavy nuclei more or less proportional to 
the SFR. Steidel et al. (1996a, 1996b) and Lowenthal et 
al. (1997) have demonstrated that galaxies with redshifts 
up to about z w 5 can be observed with optical pho- 
tometry and spectroscopy during periods of high star for- 
mation activity (resulting in high rest-frame UV emission 
which is redshifted into the optical wavelength range in the 
observers frame). Basic properties of these objects, such 
as number densities, luminosities, colors, sizes, morpholo- 
gies, star formation rates, overall chemical abundances, 
dynamics and clustering have been investigated in various 
recent papers (cf. e.g. Steidel et al. 1996b, Yee et al. 1996, 
Lowenthal et al. 1997, Pettini et al. 2000, and Leitherer 
et al. 2001). Steidel et al. (1996a, 1996b) already noted 
that the high-z galaxies show, on average, relatively weak 
metallic absorption lines, and they ascribe this finding ten- 
tatively to a lower metal content (a suggestion which has 
later been reiterated in several subsequent papers). On 
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Fig. 1. Examples of the low-resolution galaxy spectra obtained with FORS. The redshift of the objects is increasing 
from panel a) to c) . The dotted line indicates the noise level which (due to the night sky spectrum and the wavelength 
dependent instrumental efficiency) varies with wavelength. The S/N of these spectra is about < 16, 43 and 12 for 
object 7029, 5903 and 6063, respectively. The object designations correspond to the catalog of Heidt et al. (2002). For 
orientation the expected positions of selected spectral lines are indicated by vertical dashed lines. The two apparent 
emission features (marked by arrows) longward of the C iv absorption feature are artifacts and resulted from the 
increased noise at this wavelengths due to telluric OH band edges and atmospheric absorption. For the same reason 
Fig. lc becomes unreliable beyond 7000 A. 
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the other hand most of the published spectra of high-z 
galaxies are not of sufficient S/N to provide qualitative 
information on the metal content. Therefore, we obtained 
new high S/N spectra of galaxies with z < 3.5 with the 
aim of studying the chemical evolution of starburst galax- 
ies at high redshifts. 

In the present paper we describe results on the C IV 
absorption line strength and their interpretation in terms 
of chemical evolution with cosmic age at redshifts < 
z < 3.5. Most of the new spectra were obtained with the 
FORS instruments at the ESO VLT in the course of a pho- 
tometric and spectroscopic study of distant galaxies in the 
FORS Deep Field (FDF) (Appenzeller et al. 2000, Heidt 
et al. 2001, Bender et al. 2001). Due to its combination of 
depth and a (compared to the HDFs) relatively large area 
the FDF is particularly well suited for statistical studies 
of high-redshift galaxies. In the present investigation we 
restrict ourselves to redshifts z < 3.5 since at higher red- 
shifts the position of the redshifted C IV resonance lines 
tends to coincide with strong OH night sky lines. Hence, 
for an accurate sky subtraction a higher spectral resolu- 
tion or longer exposure times than we could achieve so 
far would be needed. In Sect. 2 we describe the sample 
selection and the observations. In Sect. 3 we present and 
discuss the measurement of the C iv (and Si iv) equivalent 
widths, from which we estimate the metallicities of the in- 
vestigated objects in Sect. 4. In Sect. 5 we compare our 
results with data available in the literature, in Sect. 6 we 
draw our conclusions. 

2. Sample selection and observations 

The data set used in this study combines spectra of high- 
redshift galaxies observed with the ESO VLT and spectra 
of local starburst galaxies taken from the IUE archive. 51 
of the high-redshift spectra were selected (according to the 
criteria listed below) from about 300 low-resolution spec- 
tra observed mainly during the spectroscopic observing 
runs of the FDF program (in 3 nights in Sept. and Oct. 
2000 and 3.5 nights in Oct. and Dec. 2001). For these 
observations we used FORS1&2 at the VLT in MOS and 
MXU mode with a slit width of 1" and the 1501 grism. (For 
instrumental details see the FORS Manual at the ESO 
web page www.eso.org). A few additional spectra had al- 
ready been obtained during the commissioning phases of 
FORS1&2 in 1998, 1999 and 2000 using the same spec- 
troscopical setup as described above. All spectra cover a 
spectral range (in the observer's frame) from about 3400 
A to about 10000 A with a spectral scale of 5 A/pixel 
and a spectral resolution of about 200. Although spectra 
of galaxies as faint as I = 26.0 mag have been observed 
successfully in the FDF, in the present investigation only 
galaxies with I< 24.5 mag were included. Depending on 
the objects' magnitude and the seeing conditions (Of! 7 on 
average) the integration times ranged between 2 and 12 
hours. The data reduction (bias subtraction, fiatfielding, 
cosmic ray elimination, sky subtraction, wavelength cali- 
bration, etc.), was performed using standard MIDAS rou- 



tines. A detailed description of the FDF spectroscopic ob- 
serving program and the data reduction procedures will 
be presented by Noll et al. (2002). 

With respect to the photometric redshift catalogue of 
the FDF (see Bender et al. 2001) our spectroscopic sam- 
ple is complete to about 85% for our limiting magni- 
tude (I < 24.5 mag) for photometric redshifts in the range 
2.2 < Zphot < 3.5. The distribution of the spectroscopic 
redshifts is in good agreement with the photometric red- 
shift distribution of the FDF which has peaks at redshifts 
of around 2.4 and 3.4 (see Fig. 3 in Bender et al. 2001). 
Exceptions are the redshifts in the range 1.4 < z < 2.2 
where a lack of strong spectroscopic features in our ob- 
served wavelength range makes a reliable spectroscopic 
redshift determination rather difficult, resulting in an ar- 
tificial low number of objects in our spectroscopic sample. 

For the present investigation we selected those FDF 
galaxies showing absorption line spectra with an ade- 
quate S/N (> 10 per resolution element) for a meaning- 
ful quantitative analysis of the Civ resonance doublet. 
All these galaxies show typical starburst characteristics in 
their spectra such as intense (rest frame) UV continua and 
highly ionized metal absorption lines. Three examples of 
FDF spectra are displayed in Fig. pi The 10 high-redshift 
galaxies which are dominated by their Ly Q emission were 
not included in our study. This leads to a limitation of 
the sample to galaxies with z < 4, since most of the few 
FDF galaxies with larger redshifts (up to 5.0) observed so 
far show strong Lya emission (while pure absorption line 
spectra dominate at lower redshifts; see Noll et al. 2002). 
Also excluded from our study were 4 objects that show 
absorption lines with clear emission components, forming 
P-Cygni profiles. (The apparent emission peak redwards 
from the C iv absorption feature in Fig. lc is not a P Cygni 
emission component. The absorption line is unshifted and, 
as pointed out in the caption, the apparent emission com- 
ponent is an artifact). 

In order to enlarge our sample somewhat we added 
6 additional FORS spectra (matching the criteria listed 
above) which had been observed with the same setup dur- 
ing the FORS commissioning runs (and are now available 
from the VLT archive). Four of these additional spectra 
were selected among the gravitationally amplified galaxies 
behind the cluster 1E0657-558 (Mehlert et al. 2001). We 
further included two spectra from the HDF-S and AXAF 
Deep Field follow-up studies (Cristiani et al. 2000), which 
met our criteria. 

The comparison sample of IUE low resolution spectra 
from the IUE archive^] consists of 36 local (z m 0) starburst 
galaxies investigated by Heckman et al. (1998)^. All these 
IUE spectra were obtained with the Short Wavelength 
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Prime Camera in the low dispersion mode and, therefore, 
cover a similar rest-frame spectral range with a slightly 
better spectral resolution as our FORS spectra of high 
redshift objects. The IUE spectra were reduced using the 
pipeline provided by the archive and smoothed to attain 
the same spectral resolution as the VLT spectra. 

3. Equivalent widths 

In an earlier investigation, based on smaller samples of 
high-redshift galaxies, we noticed an apparent anticorre- 
lation between redshift and the strength of the C iv 1550 A 
doublet (Mehlert et al. 2001, 2002). According to Walborn 
et al. (1995) high-excitation lines like Civ and Si IV are 
produced mostly in stellar photospheres and winds and 
their strengths depend sensitively on the stellar metallic- 
ity. Although in a few cases a non-negligible (< 50%) con- 
tamination of the Civ and Siiv features by interstellar 
absorption could not be excluded, Heckman et al. (1998) 
find for their sample of 45 nearby starburst galaxies that 
the C iv and Si iv absorption is normally produced by pho- 
tosphcric and stellar wind lines of the unresolved stars. 
We found further evidence for a close relation between 
the strength of these resonance lines and the metallicity 
of the observed starburst galaxies by measuring the equiv- 
alent widths of the C iv 1550 A and Siiv 1398 A doublets 
in synthetic spectra of starburst galaxies with different 
metallicities taken from Leitherer et al. (2001) (see Fig. ||). 
According to Fig. ||for ages > 10 Myr the measured equiv- 
alent widths depend strongly on the metallicity but are 
almost independent of the age of the starburst. Therefore, 
we measured for those galaxies with reliable spectroscopic 
rcdshifts z > 1.35 (i.e. galaxies where the C iv doublet was 
redshifted into our observed spectral range) the rest-frame 
equivalent widths Wq of this feature, which is defined by 

Here we used following approximation: 

i r X2 

Wo = AA - — — / S(X)dX with 
k(Ao) J xi 

. , AA AA 
Ai — A — ; A 2 — A H — — . 

(2) 

The central wavelength for Civ is Ao = 1549.5 A and 
for the width of the line window we chose AA = 30 A. 
The continuum flux at the central rest frame wavelength 
Ao was approximated by the mean flux within two well 
defined continuum windows, one on each side of the line 
window. Each of these continuum windows has a width 
of 75 A and is separated from the line window by 5 A. 
Since the correct continuum determination is critical for 
the resulting Wq, its level was checked interactively. In 
particular it turned out that the influence of the B-Band, 
that lies in the right continuum window of the highest-z 
objects is negligible with respect to the resulting Wq. 
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Fig. 2. Measured Civ (a) and Siiv (b) equivalent width 
of the synthetic spectra of Leitherer et al. (2001) as a func- 
tion of the starburst age. The model spectra are based on 
continuous star formation (1 M Q /yr) assuming the model 
parameter M up = 100 M and aiMF = 2.35. Circles and 
triangles correspond to solar and LMC metallicity, respec- 
tively. 



For objects with z > 1.7, where in addition the Siiv 
doublet (Aq = 1398.3 A) became visible, we also measured 
Wo (Siiv) with the same bandwidth and continuum win- 
dow definitions as for C iv. The same measurements were 
also carried out in the IUE spectra of the comparison sam- 
ple of z w starburst galaxies. Statistical mean errors of 
the individual Wq measurements were calculated from the 
S/N of the individual spectra and the errors of the con- 
tinuum fits (< 10%). Equivalent width measurements in 
the spectra of faint objects can be affected severely by 
errors in the sky background subtraction. Therefore, dur- 
ing the reduction of our FORS spectra a particular effort 
was made to keep these errors low. Various tests showed 
that for all FORS spectra used in this study the errors in 
determining the continuum level remained below 5% (see 
Noll et al. 2002). Hence, systematic errors in the measured 
equivalent widths due to an incorrect sky subtraction are 
well below the statistical errors in most cases. 

Another source of systematic errors in our Wq mea- 
surements is our low spectral resolution, which in most 
cases did not allow us to resolve the line profiles. However, 
since equivalent widths measurements of strong isolated 
lines are in principle independent of the spectral resolu- 
tion and since all our conclusions are based on differences 
between measurements carried out with the same proce- 
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Table 2. Averages and mean errors of the measured C iv A 
1550 rest frame equivalent widths within selected redshift 
bins for all galaxies listed in Table |l|. N is the number of 
objects within each bin. 
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Fig. 3. Measured C iv A 1550 rest-frame equivalent widths 
as a function of redshift. The various symbols have the fol- 
lowing meaning: Open star at z w 0: Average and lcrrms 
scatter for the 36 (Heckman et al. 1998) local starburst 
galaxies. Open triangles: FDF galaxies. Filled triangles: 
Galaxies in the field of the cluster 1E0657-558 (Mehlert et 
al. 2001). Filled circle and square: Galaxies in the HDF-S 
and the AXAF Deep Field, respectively (Cristiani et al. 
2000). 
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Fig. 4. Averages of the measured C iv A 1550 rest 
frame equivalent widths within selected redshift bins (see 
Table ||) as a function of redshift for all galaxies shown in 
Fig. H (asterisks) . The bars indicate the mean errors of the 
averages. Open triangles show all FDF galaxies brighter 
than Mb = —22.28 mag, open squares all FDF galaxies 
with -21.52 mag < M B < -20.38 mag (see Sect. § for 
discussion) . 



dure in spectra of the same resolution, these systematic 
errors are expected to cancel out and therefore are not 
expected to affect the results of this paper significantly. 
On the other hand, since the C iv doublet is not a truly 
isolated feature, our numerical results for the equivalent 
widths should not be directly compared to results obtained 
from spectra with a different spectral resolution. 

The results of our Wo(C iv) measurements are listed in 
Table [j] and plotted in Fig. ||. In order to avoid a crowding 
of data points at z ss 0, we plotted for the local (IUE) 
starburst galaxies only the average value and indicate the 
lcr scatter of the individual values by a bar. For the high- 
redshift galaxies the individual data points and their mean 
errors are given. 



z interval 


< z > 


< Civ > 

[A] 


m.e.(C iv) 

[A] 


N 


0.00 - 0.02 


0.01 


6.31 


0.41 


36 


1.00 - 1.99 


1.52 


4.37 


0.41 


6 


2.00 - 2.49 


2.32 


4.51 


0.27 


32 


2.50 - 2.99 


2.71 


3.77 


0.53 


7 


> 3.00 


3.28 


2.33 


0.22 


12 



As demonstrated by Fig. || our high-redshift galaxies 
with z < 2.5 show about the same average C IV equivalent 
widths and about the same scatter around the average as 
the local starburst galaxies. However, for redshifts larger 
than about 2.5 the average Civ equivalent widths and 
their scatter clearly decrease with z in our sample. As 
described in Sect. || this decrease is not driven by any se- 
lection effect since the spectroscopic redshift distribution 
of the included FDF galaxies is in good agreement with 
the photometric redshift distribution. 

Fig. H obviously confirms the effect suspected by 
Steidel et al. (1996a) quantitatively. In order to estimate 
(in view of the observed scatter) the statistical significance 
of the effect, we calculated averages and their mean er- 
rors of the Wo(Civ) values for selected redshift bins. The 
results are listed in Table ^ and plotted in Fig. ^. The 
table confirms that there is no statistically significant dif- 
ference between the results for the first three bins while 
the difference between the local sample and our starburst 
galaxies with z > 3.0 is highly significant (> 9er). Tests 
with other bin sizes and binning intervals showed that the 
high significance of the result persists for any reasonable 
bin distribution. 

In Fig. H we present analogously to Fig. || the observed 
Si iv equivalent width values as a function of redshift. 
Similarly as in the case of the C iv doublet the average 
Si iv strength does not change with redshift for z < 2.5. 
However, unlike the average C IV strength the average Si iv 
Wo values remain at the local value even beyond z > 2.5. 
As a result, the ratio between the Siiv and C iv resonance 
doublets, which is practically constant for low z, varies for 
high redshifts in our sample. This is demonstrated quan- 
titatively by Table | and Fig. g. 

As noted above, for local starburst galaxies the C iv 
doublet as well as the Siiv doublet are good metallicity 
indicators. Since both elements are produced during the 
evolution and explosion of massive stars, a greatly differ- 
ent relative chemical abundance of C and Si in the high-z 
starburst galaxies appears very unlikely. In spite of the un- 
satisfactory state of present SN II models and the remain- 
ing large uncertainties concerning the intermediate mass 
element yields for different initial stellar masses, it seems 
very difficult to enhance Si production relative to C. On 
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Table 1. Measured Civ and Siiv rest-frame equivalent widths for the 57 high-z galaxies. The absolute B-magnitudes 
derived and discussed in Sect. H are also listed. 



No. 


z 


I 


Wo (Civ) 


d W (Civ) 


Wo (Siiv) 


dWo(Siiv) 


M B 






[mag] 


[A] 


[A] 


[A] 


[A] 


[mag] 


FDF-1208 


2.18 


23.68 


3.92 


0.41 


3.31 


0.44 


-22.45 


FDF-1331 


3.39 


23.89 


1.80 


0.90 


_ 


_ 


-23.30 


FDF-1555 


3.26 


23.88 


1.16 


0.74 


1.29 


0.69 


-22.36 


FDF-1578 


2.71 


24.25 


3.545 


0.90 


2.97 


0.81 


-21.98 


FDF-1691 


2.34 


23.89 


3.65 


0.57 


4.30 


0.60 


-21.53 


FDF-1709 


1.67 


24.33 


6.07 


0.75 


— 


— 


-20.46 


FDF-1744 


2.37 


24.10 


2.34 


0.58 


4.27 


0.56 


-22.21 


FDF-1922 


1.83 


23.36 


3.62 


0.28 


0.95 


0.38 


-21.64 


FDF-2033 


2.75 


24.08 


5.51 


0.65 


4.05 


0.52 


-21.62 


FDF-2274 


2.25 


23.34 


1.62 


0.31 


2.66 


0.33 


-21.66 


FDF-2418 


2.33 


23.16 


6.57 


0.47 


6.36 


0.46 


-23.19 


FDF-2495 


2.45 


23.31 


3.98 


0.32 


3.24 


0.32 


-22.32 


FDF-2636 


2.25 


23.43 


5.21 


0.63 


6.41 


0.70 


-22.77 


FDF-3005 


2.25 


23.51 


6.95 


0.38 


8.18 


0.43 


-22.63 


FDF-3163 


2.44 


23.35 


4.80 


0.33 


5.45 


0.33 


-22.97 


FDF-3173 


3.27 


23.91 


2.59 


0.48 


4.24 


0.44 


-22.51 


FDF-3300 


2.37 


23.91 


2.14 


0.42 


2.16 


0.41 


-21.79 


FDF-3374 


2.38 


23.34 


5.05 


0.30 


5.21 


0.27 


-22.65 


FDF-3810 


2.37 


22.67 


4.95 


0.25 


5.59 


0.26 


-23.18 


FDF-3874 


2.48 


23.30 


2.66 


0.43 


3.83 


0.43 


-23.15 


FDF-3875 


2.24 


24.53 


4.19 


0.51 


3.38 


0.52 


-20.73 


FDF-3958 


2.13 


23.87 


4.40 


0.53 


1.43 


0.57 


-20.98 


FDF-3999 


3.39 


24.00 


3.74 


0.65 


7.89 


0.56 


-22.70 


FDF-4049 


1.48 


23.00 


4.42 


0.53 


— 


— 


-21.76 


FDF-4795 


2.16 


23.31 


5.80 


0.36 


5.93 


0.39 


-22.35 


FDF-4871 


2.47 


23.39 


7.14 


0.35 


6.05 


0.34 


-22.62 


FDF-4996 


2.03 


23.25 


3.35 


0.37 


1.19 


0.44 


-21.77 


FDF-5058 


2.03 


23.34 


4.82 


0.25 


3.40 


0.27 


-21.57 


FDF-5072 


1.39 


22.45 


3.27 


0.57 


_ 


_ 


-21.88 


FDF-5135 


2.34 


23.62 


2.31 


0.71 


1.99 


0.77 


-22.73 


FDF-5152 


1.37 


22.65 


3.91 


0.50 


— 


— 


-21.55 


FDF-5165 


2.35 


23.26 


5.73 


0.55 


6.02 


0.53 


-23.02 


FDF-5190 


2.35 


24.39 


2.73 


0.68 


2.65 


0.64 


-22.74 


FDF-5215 


3.15 


22.98 


2.54 


0.47 


2.45 


0.40 


-23.18 


FDF-5227 


2.40 


23.85 


4.73 


0.72 


2.02 


0.79 


-21.79 


FDF-5504 


3.38 


23.63 


3.29 


1.00 


7.19 


0.73 


-23.65 


FDF-5550 


3.38 


23.12 


2.44 


0.41 


5.23 


0.31 


-23.23 


FDF-5903 


2.77 


22.33 


4.02 


0.21 


4.32 


0.16 


-23.23 


FDF-6024 


2.37 


22.00 


4.93 


0.20 


5.71 


0.19 


-23.34 


FDF-6063 


3.40 


22.56 


1.36 


0.49 


4.24 


0.42 


-23.35 


FDF-6069 


2.68 


24.22 


5.36 


0.81 


3.45 


0.70 


-21.74 


FDF-6287 


2.68 


24.11 


1.61 


0.82 


- 


— 


-22.04 


FDF-6372 


2.35 


23.38 


3.20 


0.37 


3.43 


0.35 


-22.18 


FDF-6407 


2.16 


23.59 


5.29 


0.59 


4.26 


0.72 


-22.31 


FDF-6864 


1.39 


23.41 


4.93 


0.71 


- 


- 


-20.87 


FDF-6934 


2.44 


22.90 


4.71 


0.56 


3.01 


0.60 


-23.05 


FDF-6947 


2.36 


23.83 


5.91 


0.42 


4.90 


0.45 


-21.94 


FDF-7029 


2.37 


23.63 


6.88 


0.42 


5.82 


0.42 


-23.07 


FDF-7307 


2.44 


24.07 


2.51 


0.60 


3.05 


0.62 


-21.20 


FDF-7342 


2.37 


23.80 


6.13 


0.76 


4.58 


0.75 


-22.06 


FDF-7539 


3.29 


23.51 


1.92 


0.46 


4.22 


0.41 


-22.41 


ES0657-A 


2.34 


24.50 


5.88 


0.48 


3.98 


0.45 




ES0657-C 


3.08 


24.89 


2.30 


1.00 


7.07 


0.62 




ES0657-J 


2.61 


22.98 


3.81 


0.42 


2.41 


0.26 




ES0657-Core 


3.24 


24.31 


2.01 


0.63 


3.45 


0.53 




HDFS-047 


2.79 




2.50 


0.30 


4.40 


0.17 




AXAF-028 


3.13 




2.84 


1.00 


3.45 


0.97 
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Redshift 



Fig. 5. Measured Siiv A 1400 rest-frame equivalent width 
as a function of redshift. Symbols as is Fig. [|. 



the other hand, it is well known that individual hot stars 
of the same metallicity show a wide range of C iv to Si iv 
line ratios. Strong C IV absorption is well known to be uni- 
versally present in O stars of all luminosity classes, while 
Wo(Siiv) is luminosity dependent and decreases rapidly 
from supergiant stars to dwarfs (Walborn & Panek 1984, 
Pauldrach et al. 1990, Leitherer et al. 1995). This line 
is, therefore, used as a luminosity indicator in UV stellar 
classification schemes. Moreover, Si IV has a pronounced 
maximum in early B stars while Civ changes monotoni- 
cally with temperature. Finally the Si iv strength is more 
strongly affected by population differences (i.e. stellar age 
differences) than the Civ doublet. A scatter in these pop- 
ulation differences can easily mask any metallicity depen- 
dence in the Siiv line strength. 

Hence, assuming that the observed absorption lines are 
dominated by contributions of the stellar photospheres 
and winds, the Wo (Si iv)/Wo(C iv) ratio can, in princi- 
ple, be used to derive informations on the star formation 
history (instantaneous or continuous) and/or the stellar 
mass distribution (see also Mas-Hesse & Kunth 1991). 
Therefore, Fig. ^ can possibly be understood by assum- 
ing that at the epochs corresponding to z > 3 (i.e. dur- 
ing the first two Gyrs of the universe) instantaneous star 
bursts played an important role while "continuous star for- 
mation" is the normal mode for local and lower redshift 
starburst galaxies. However, since the different parameters 
which determine the value of Wo (Si iv)/Wo(C iv) cannot 
be disentangled reliably without including additional lines 
in the analysis, spectra of higher resolution and higher S/N 
are required to finally settle this issue. 

The investigation of the purely interstellar lines of 
lower ionisation like e.g. Sill A 1260, Ol/Sin A 1303 and 
C II A 1335 is in progress and the results will be published 
in a separate paper. 

4. Metallicities 

Since, in contrast to the variations of the Siiv line, differ- 
ences of the C iv line strength cannot be easily explained 
by population differences in the starburst galaxies, the 
observed decrease of the C IV equivalent width values for 



Redshift 

Fig. 6. Averages of the ratios between the measured Si iv 
and C IV rest frame equivalent widths as a function of 
redshift for all galaxies shown in Fig. The mean error of 
each average, calculated from the scatter of the individual 
values, is indicated by a bar. 

Table 3. Average and mean error of the ratios between 
the measured Si IV and C IV rest frame equivalent widths 
within selected redshift bins for all high-z galaxies for 
which Wo (Siiv) has been measured. N is the number of 
galaxies within each bin. 



< z > 


< Siiv/Civ > 


m.e.(Si iv/C iv) 


N 


0.00 


0.76 


0.04 


36 


2.30 


0.94 


0.06 


33 


2.72 


0.95 


0.18 


6 


3.27 


1.95 


0.22 


11 



z > 2.5 in our sample can at present only be interpreted as 
a metallicity effect. Hence, the decrease of Wo (Civ) with 
z is expected to contain information on the evolution of 
the metal content of starburst galaxies with cosmic age. In 
order to derive a more quantitative measure of the metal- 
licity evolution apparently observed in Fig. ||, we made 
an attempt to calibrate the observed W (Civ) values in 
terms of the O/H ratios. For this purpose we used the 
oxygen abundances listed in Heckman et al. (1998) for 
all the local starburst galaxies of this sample and derived 
metallicities using the relation \ogZ = 12 + \og{0/H). 

In Fig. |?] we plotted for the local starburst galaxies 
the metallicity relative to the solar value (logZo = 8.93) 
as a function of our measured Wo (Civ) values. We also 
included the theoretical values determined for synthetic 
starburst galaxies with solar and LMC metallicities taken 
from Leitherer et al. (2001). Although the scatter is rather 
large (RMS = 0.27), the plot indicates a dependence of 
log Z on Wo(C iv) which can be approximated by a linear 
relation. The best linear least square fit to these data gives 

logZ/Z = 0.13(±0.02) x Wo(Civ) - 1.1Q(±0.12). (3) 

This calibration of the C iv strength in terms of metal- 
licity should be a reasonable approximation for statisti- 
cal applications at least for local starburst galaxies and 
the objects with z < 2.5. However, in view of the popu- 
lation differences evident from the different Si iv to C IV 
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Fig. 7. Metallicity log Z/Zq in terms of oxygen abundance 
of 36 local starburst galaxies (from Heckman et al. 1998) 
as a function of the measured C iv A 1550 equivalent width 
(open stars). The typical error for \ogZ/Z Q is < 0.1 dex. 
The solid line gives the best linear fit to the data. For 
comparison the W^Civ) of synthetic starburst spectra 
(from Leitherer et al. 2001; see Fig. ||) at an age of 100 
Myr for solar (filled circle) and LMC (0.25 solar, filled 
triangle) metallicity have also been included in the figure. 




Fig. 8. Measured Civ A 1550 rest-frame equivalent 
widths of the local starburst galaxies (open stars) and 
the high-z FDF (filled symbols) versus the absolute B- 
magnitude. For the latter ones we indicated objects within 
the different redshift bins l<z<2,2<z<3 and 
z > 3 by filled circles, triangles and squares, respectively. 
The solid vertical line indicates Mb = —22.36 mag , while 
the dashed vertical lines indicate Mb = —21.52 mag and 
M B = -20.38 mag 



ratio, its applicability to the z > 2.5 objects is less clear. 
Nevertheless, because of the absence of other more reliable 
calibration procedures, and since (in view of the physics 
of hot stars) the relation between the C iv strength and 
the metallicity should not be much affected by popula- 
tion details of starbursts we will assume for the follow- 
ing that the correlation between the oxygen abundances 
and W(C iv) observed for low redshift starburst galax- 
ies is also valid at high redshifts. With this assumption 
we convert our observed C IV equivalent width values to 
metallicities using Eq. ||. In this way we obtain for our 
starburst galaxies with z > 3 (< z >= 3.24) an average 
metallicity of about 0.16 Zq and for < z >= 2.34 a value 
of 0.42 Zq. The corresponding local (z — 0) value would 
be 0.56 Zq. In terms of cosmic time scales (for a universe 
with n A = 0.7, n M = 0.3, H = 67 km sec" 1 Mpc" 1 used 
throughout the paper) this would correspond to an in- 
crease of the mean metallicity in starburst galaxies by a 
factor of 2.5 within w 1 Gyr between cosmic ages of about 
1.9 Gyrs and 2.9 Gyrs. For later epochs the data suggest 
only little further enrichment. Because of the approxima- 
tive nature of Eq. || these numbers are rough estimates 
only. Still, they agree surprisingly well with earlier theo- 
retical predictions of the cosmic chemical enrichment his- 
tory of the universe by e.g. Fritze-von Alvensleben (1998) 
and by Renzini (1998, 2000), who predicts that the metal- 
licity had been w 0.1 Zq at z = 3 and has increased to a 
value of ~ 1/3^0 in the local universe. 

Two effects may affect the validity of Eq. || at high 
redshifts: First, in the IUE spectra of the local starburst 
galaxies contributions of the Milky Way halo components 
are present, while they are absent in the high-z spectra. 
Savage & Massa (1987) showed that the Civ equivalent 
width of distant halo stars are normally < 0.5 A. This is 



one order of magnitude lower than the values measured 
in the local starburst galaxies. Hence its effect on Equ. 
3 should be on a 10% level, at most. Secondly, the spec- 
trum of the magnified high-z object MS1512-cb58 (z = 
2.727) published by Pettini et al. (2000) suggests that the 
contribution of the interstellar line to the C iv absorption 
feature increases with redshift. If this holds for all galaxies 
at similar redshifts, this could result in a general difference 
of the mean C iv equivalent width between local starburst 
galaxies and high-redshift objects, but it could not ex- 
plain the observed evolution of the C IV equivalent width 
between z = 2.5 and z = 3.5. Furthermore, by applying 
the local relation to high redshifts we would overestimate 
the mean metallicity at young epochs. 

5. Luminosity effects 

When interpreting the data described above one has to 
keep in mind that for local (z « 0) galaxies the metal- 
licities are known to depend on the galaxies' blue and 
infrared luminosities, with luminous galaxies tending to 
have higher metallicities (see e.g. Kobulnicky & Zaritsky 
1999 and Heckman et al. 1998). Since at high redshifts 
we observe only very luminous galaxies. Therefore, if a 
metallicity-luminosity correlation exists these bright ob- 
jects should be metal rich and we should find an opposite 
correlation between metallicity and redshift than the one 
detected in this work. To test whether at high redshifts a 
metallicity-luminosity correlation does exist and may af- 
fect our detected metallicity evolution we plotted in Fig. || 
the absolute B-magnitudes Mb of all high-redshift FDF 
galaxies as well as Mb for the local starburst galaxies. 
For the local objects the Mb was taken from Heckman et 
al. 1998 (transformed to the cosmology used in this pa- 
per), for the FDF galaxies we computed Mb as follows: 
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We derived the best fitting SED, scaled to the total 1 flux 
derived by SExtractor (FLUX.AUTO) as determined by 
our photometric redshift code (see Bender et al. 2001). 
Then this SED was transformed to z = (using the ob- 
served spectroscopic redshift) to derive the rest-frame B- 
magnitude of the galaxies. Since for the redshift range in 
question the measured J and K bands bracket the rest- 
frame B, this procedure is nearly equivalent to an inter- 
polation, minimizing the uncertainties in the K correc- 
tions. A detailed description of the method can be found 
in Gabasch et al. (2002). Using the photometric instead of 
the spectroscopic redshifts would produce a typical vari- 
ations of ±0.2 mag. Absolute magnitudes were derived 
assuming the cosmology parameter Hq — 67, fi m = 0.3, 
17a = 0.7. Our Mb have been corrected for foreground 
Galactic extinction but not for airy internal extinction in 
the starburst galaxies. From Fig.|| we see that the local 
starburst galaxies indeed show the expected correlation 
between Wo(Civ) and the luminosity. On the other hand, 
for the high-redshift galaxies we cannot determine whether 
a metallicity-luminosity relation does exist or not, since we 
do not have any faint objects in our high-z sample. But 
it is evident that the high-redshift galaxies are on aver- 
age overluminous for their metallicities compared with lo- 
cal starburst galaxies. This agrees well with earlier results 
from Pcttini et al. (2001) and Kobulnicky & Koo (2000) 
who find this trend for Lyman-break galaxies. Hence, if a 
metallicity-luminosity relation does exist at high redshifts, 
our data suggest that it has a clear offset to the local cor- 
relation, which seems to evolve with redshift. Moreover, 
from FigJ| it is obvious that for the high-redshift galax- 
ies there is no correlation between the measured Wo(C IV) 
and the luminosity that could cause the correlation with 
z found in this paper. 

Since at high redshifts we do not have any faint ob- 
jects in our sample, while in the local universe we do not 
find bright starburst galaxies, we have to make sure that 
our detected metallicity evolution with redshift is not pro- 
duced by comparing different objects at different redshifts. 
For that reason we separately investigated all galaxies, 
which are brighter than the faintest one at z > 3 (which 
is Mb = —22.36 mag; solid line in Fig. |J). In our sample we 
only find galaxies brighter than this limit for z > 2. Their 
average values of the measured Wo (Civ) and the mean 
error at redshift 2.4 and 3.3 as well as the single Wo(Clv) 
measurement for this brightest local galaxy are addition- 
ally indicated by open triangles in FigJJ. Obviously these 
subsample show the same trend with decreasing redshift as 
the total galaxy sample at z > 2. Furthermore we investi- 
gated all galaxies fainter than Mb = —21.52 mag (bright- 
est local galaxy) and brighter than Mb = —20.38 mag 
(faintest galaxy with z > 1). The average values of the 
measured Wo(Civ) and the mean error at redshift 0, 1.5 
and 2.4 are also indicated in Fig.f| by open squares and 
again show the same trend with decreasing redshift. From 
these test we conclude that the observed dependence of 
Wo (Civ) on redshift is not caused by a luminosity effect. 
Moreover the two open symbols in Fig. ^ for z rj 2.4 in- 




o - 



40 40.5 41 

log(L uv [ergs- 1 A-']) 

Fig. 9. Measured Civ A 1550 rest-frame equivalent 
widths of the high- 2 FDF galaxies as a function of the UV 
luminosity Ltjv, derived from the flux between 1432 A and 
1532 A (as defined by Kinney et al. 1993) 

dicate that a metallicity-luminosity correlation also exists 
at this redshift. 

The following additional selection effects could be 
present (and possibly weaken) the observed correlation be- 
tween metallicity and redshift at high-z: ft could, in princi- 
ple, be possible that at high-z we preferentially see objects 
with low internal extinction, which have low dust content 
and hence low metallicity. In this case we would expect 
to find a negative correlation between the UV luminosity 
of our galaxies and their metallicity. To test whether this 
correlation is present in our high-z galaxies we calculated 
the UV luminosity as follows: 

log(L uv ) = log < F >«g 

+ log(4 7r(10pc) 2 ) 

(4) 

where <F > 1432 is the mean flux between 1432 A and 
1532 A (as defined by Kinney et al. 1993), (m-M) is the 
distance modulus calculated using the cosmological pa- 
rameters mentioned above. We did not correct for intrin- 
sic reddening since such a correction would have involved 
considerable uncertainties and is not needed for the the 
present test. From Fig. [|, where we plotted the C IV equiv- 
alent widths for the high-z FDF galaxies as a function of 
LuV; w e see that no such correlation is evident for our 
high-z objects. Hence our galaxy sample seems not to be 
affected by this selection effect. 

6. Comparison with literature data 

Although most published spectra of high redshift galaxies 
lack the S/N required to carry out a study of the type 
presented here we tried to compare our results with the 
limited information available on this subject in the liter- 
ature. For this purpose we measured the Civ equivalent 
width for all high-z galaxy spectra published in the pa- 
pers listed in the caption of Table ^. As in our samples, 
spectra with strong Lya emission were disregarded, and 
only equivalent width values with estimated (rest frame) 
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Fig. 10. Comparison of the observed ClV A 1550 rest 
frame equivalent widths as a function of redshift for 
our galaxy sample (asterisks) with the mean values of 
Wo (Civ) derived from high-z galaxy spectra from Pettini 
et al. 1998 (P98; open circle), Steidel et al. 1996a (S96; 
open square), Steidel et al. 98 (S98; open hexagons) and 
Lowenthal et al. (L98; open diamond). The bars denote 
mean errors. Single measurements from galaxies observed 
by Pettini et al. 2000 (POO; filled circle), Yee et al. 1996 
(Y96; filled square) and Trager et al. 1997 (T97; filled tri- 
angle) are also shown. Note that the Y96 and POO observed 
the same object namely the well known lensed galaxy 
MS1512-cb58. 

Table 4. Line 1-4: Average of the measured C iv A 
1550 rest-frame equivalent widths for high-z galaxies ob- 
served by Pettini et al. (1998; P98), Lowenthal et al. (1997; 
L97) and Steidel et al. (1996a & 1998; S96 & S98). Line 
5-7: Measured C iv A 1550 rest frame equivalent width 
for MS1512-cb58 observed by Yee et al. (1996; Y96) and 
Pettini at al. (2000; POO) and for object DG-433 observed 
by Trager et al. (1997; T97). A * indicates that the equiva- 
lent widths were measured using the same method applied 
to the the high-z galaxies presented in this work. A f in- 
dicates that the equivalent widths were measured from 
enlarged tracings. 



< z > 


< Civ > 

[A] 


m.e.(C iv) 

[A] 


N 


Reference 


2.92 


3.55 


1.23 


5 


P98* 


3.05 


1.80 


0.39 


5 


L97 f 


3.09 


2.28 


0.59 


4 


S98 f 


3.22 


2.20 


0.50 


2 


S96 1 " 


< z > 


Wo(Crv) 

[A] 


dWo(Civ) ) 

[A] 


N 


Reference 


2.73 


5.14 


1.20 


1 


POO* 


2.73 


2.91 


1.11 


1 


Y96* 


3.34 


2.30 


0.70 


1 


T97 f 



mean errors < 1.0 A were used for the comparison. The 
number of the (for our purpose) usable spectra of each 
publication is listed in column 4 of Table |[ (Tests showed 
that including less accurate data lead to similar results, 
but with much larger statistical errors). For the objects 
observed by Yee et al. (1996) and by Pettini et al. (1998, 
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Fig. 11. Mean of the estimated metallicities as a func- 
tion of redshift for all galaxies shown in Fig. ||. The oxy- 
gen abundances and the statistical uncertainties derived 
from i?23 ratio for 5 Lyman-break galaxies investigated 
by Pettini et al. (2001) and Teplitz et al. (2000) are indi- 
cated by the vertical lines: solid line: upper branch results; 
dashed line: lower branch result. Solar and LMC metallic- 
ities are indicated by the horizontal dotted lines. 



2000) the author kindly made their spectra available to us 
in electronic form. Hence, to measure Wo(Civ) in these 
literature objects we were able to apply exactly the same 
procedure as used for the galaxies presented in this work. 
For the objects investigated by Steidel et al. (1996a, 1998), 
Lowenthal et al. (1997) and Trager et al. (1997) we mea- 
sured the Wo (Civ) from enlarged tracings. We tested the 
reliability of measuring Wo(C iv) from tracings using some 
of our own high-z galaxies. The difference between the two 
measurement methods turned out to be < 10%. Although 
the spectral resolutions of the different investigations are 
not exactly the same, the individual resolutions are suffi- 
ciently close to allow a direct comparison within the accu- 
racy needed here. Tableland Fig. [To|show that the aver- 
age C iv values derived from the literature high-z spectra 
are in reasonably good agreement with the mean values 
derived for the galaxies investigated in this work. However, 
the scatter of the literature data at high redshift is larger. 
Although this larger scatter is presumably dominated by 
the on average low S/N of the literature spectra, we can- 
not exclude that environmental effects may influence the 
evolution of the C iv strength at high redshifts. 

The best investigated individual high-z galaxy is, so 
far, the gravitationally magnified object MS1512-cb58 
(z = 2.727) (cf. e.g. Yee et al. 1996, Seitz et al. 1998, 
Pettini et al. 2000, Teplitz et al. 2001, Savaglio et al 2002). 
By measuring the C iv equivalent width on low resolution 
spectra and using Eq. |^ we obtain for this galaxy a metal- 
licity of O.4Z and 0.2Z Q from POO's and Y96's data, re- 
spectively. Within our error limits these values are in good 
agreement with the result of Pettini et al. (2000) (who 
derive 0.25Z Q by comparing the galaxy spectrum with 
synthetic starburst galaxy spectra from Leitherer et al. 
(2001)) and Teplitz et al. 2001 (who found O.32Z from 
measuring its oxygen abundance on N1R spectra using 
the strong line index R23 which relates (O /H) to the rela- 
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tive abundance of [Oil], [OIII] and H/3). This comparison 
seems to support our assumption that our calibration of 
the C iv strength in terms of metallicity is applicable to 
high redshift objects at z = 2.7, although the redshift of 
MS1512-cb58 is too small to estimate the accuracy of the 
method for the interesting z > 3 objects. 

For four further high-z Lyman-break galaxies Pettini 
et al. (2001) determined the oxygen abundance from NIR 
spectra using the strong line index -R23. Since the relation 
between R23 and O/H has an upper and a lower branch 
these abundances show the well-known two-value ambi- 
guity. Hence these results do not provide a reliable test 
of our conclusions. Nevertheless, in Fig. [ll] we plot the 
allowed ranges of oxygen abundance for these 4 Lyman- 
break galaxies as well as for MS 1512-cb58 together with 
the metallicities of our starburst galaxies (as derived from 
the C iv strength via the calibration described above) as a 
function of redshift. This comparison shows that all data 
are at least mutually compatible, although for the two 
highest- z galaxies from Pettini et al. (2001) only the lower- 
branch results are in reasonable agreement with a strong 
increase of metallicity from redshift w 3.2 to w 2.3 sug- 
gested by our results. 

Our results are also in line to those obtained by de 
Breuck et al. (2000), who find a qualitative increase of 
metallicity from higher to lower redshift for a sample of 
high-z radio galaxies. Furthermore Pettini et al. (1997) 
and Savaglio et al. (2000) report on evidence for a grad- 
ual chemical enrichment of the gas producing the damped 
Lycv lines in QSO spectra, although their trends are only 
weakly significant. Compared to Savaglio et al. (2000) we 
find a zero point offset of the metallicity-redshift relation 
of about 0.7 in A log Z at z — 2.5. Such a difference is not 
unexpected since the metal absorbers in damped Lya sys- 
tems most likely sample the outermost regions of galaxies 
and therefore a different environment than the dense inter- 
stellar matter of which the massive stars seen in starbursts 
have been formed. 

7. Conclusions 

Our study shows that the FDF starburst galaxies at 
z ~ 3.2 have on average significantly lower Civ equiva- 
lent widths than starburst galaxies at lower redshifts. In 
view of the known close relation of the C IV strength to 
the metallicity in local starburst galaxies, it appears likely 
that this effect is due to a significant evolution of the av- 
erage metallicity in such objects at high redshifts. Using 
data from well studied local starburst galaxies we cali- 
brated the C iv strength in terms of the heavy element 
content of these objects. Assuming that this calibration 
is applicable to high redshift starburst galaxies we find 
that the mean cosmic metallicity as observed in starburst 
galaxies has increased significantly between the cosmic 
epochs corresponding to z « 3.2 and « 2.3 (when the 
universe was between 2 and 3 Gyrs old). If this interpre- 
tation of the increase of the C IV absorption strength with 
decreasing redshift is correct, an intense phase of star for- 



mation and evolution of massive stars must have occurred 
during this period. At lower redshifts (z < 2.5) our data 
indicate little further increase of the average metallicity of 
starburst galaxies. Hence, the further cosmic chemical en- 
richment seems to have been insignificant during the last 
11 Gyrs. The metallicity evolution indicated by our data is 
in reasonable agreement with published theoretical chem- 
ical enrichment models due to star formation at early cos- 
mological epochs. Our results are also in agreement with 
a metallicity evolution found in high-z radio galaxies (de 
Breuck et al. 2000) and the tentative evidence for a grad- 
ual chemical enrichment of the gas producing the damped 
Lya lines (Pettini et al. 1997, Savaglio et al. 2000). 

Since the FDF observations provide a fairly complete 
sample of the bright starburst galaxies in the observed 
direction and redshift ranges, the observed chemical evo- 
lution should be characteristic of the cosmic volumes 
with the most intense star formation at the correspond- 
ing epochs. However, it is also clear, that our results do 
not apply to all objects and volumes at a certain red- 
shift or epoch. It is, e.g., well known that the BLR gas of 
high-redshift QSOs is characterized by high metallicities 
and that no significant chemical evolution is observable in 
these objects up to at least z = 5 (see e.g. Hamann and 
Ferland, 1999, Dietrich et al. 1999, Dietrich and Wilhelm- 
Erkens, 2000). Obviously, in the environment of these 
early QSOs much star formation and stellar evolution 
must have taken place at epochs corresponding to even 
higher redshifts. It is also known, however, (and confirmed 
from the FDF survey) that high-z starburst galaxies are 
much more frequent than bright high-z QSOs. Hence, the 
starburst galaxies are expected to be a more representa- 
tive tracer of the history of the overall cosmic chemical 
evolution than the QSOs. 

From our data we cannot determine whether at high 
redshift a metallicity-luminosity relation does exist, since 
we do not have any faint objects in our high-z sample. But 
it is evident that the high-redshift galaxies are on average 
overluminous for their metallicities compared with local 
starburst galaxies. This trend is also found by Pettini et 
al. (2001) and Kobulnicky & Koo (2000) for Lyman break 
galaxies. From tests on various subsamples we find that 
the observed dependence of Wo (Civ) on redshift is not 
caused by a luminosity effect. Furthermore we showed that 
this dependence is also not cause by the possible selection 
effect of preferentially observing galaxies with low dust 
content and hence low metallicity at high redshifts. 

Differences in the Si iv to C iv ratios between local 
and high-z galaxies in our sample suggest differences in 
the population and star formation history in the galaxies 
with z > 3. Short bursts of star formation may have been 
more important (relative to periods of "continuous star 
formation" ) at these early epochs. 
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